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I. INTRODUCTION 

In a recent paper |TJ a statistical theory of the bound states of atomic nuclei has been 
developed, where a new definition of the nuclear mean field (MF) and of the fluctuations 
around it (the "error" ) are given. 

The aim of this theory is to provide an alternative route to overcome the problem pre- 
sented by the strong forces acting among the nucleons inside nuclei at short distances. These 
indeed entail very large matrix elements between the ground and the highly excited nuclear 
states, thus mixing into the former quite complex configurations. For this physics, in spite 
of the many approaches worked out in the past, perturbative treatments are of little value. 

Hence, the MF is obtained in [0] through an energy average that smoothes out short time 
events. Furthermore, the matrix elements of the residual interaction (the one remaining after 
the average interaction is removed) are assumed to be random with vanishing average value. 
These two elements are sufficient to obtain the MF and the "error" , which of course should 
be small for the MF to be meaningful. 

Concerning the error, while its average is zero, the average of its square is not and, in 
our framework, it can be computed. Indeed, we have organized the energy fluctuations 
in classes of excited states of increasing complexity, each class yielding a contribution to 
the error. Notably, a criterion for the convergence of the expansion of the latter in the 
complexity of the states has been obtained. 

The above outlined statistical framework has been applied in |T]] to the problem of nuclear 
matter, in order to provide an appreciation of the results the approach leads to. In carrying 
out this analysis, however, rough approximations were made, which might cast some doubts 
on the validity of our findings. Prominent among the latter is the one related to the dramatic 
reduction of the strength of the residual interaction, with respect to the bare one. Also the 
very fact that a statistical computation of the ground state energy of nuclear matter yields 
a reasonable result with a spectroscopic factor varying between consistent values across the 
energy range of the fluctuation is relevant. 

In this paper we remove the approximations referred to above and explore whether our 
previous findings survive the test of a more accurate analysis. At the same time we examine 
the dependence upon the two inputs required by our scheme, namely the Hartree-Fock 
(HF) energy and the average energy of the two particle-two hole (2p-2h) excitations. Of 
significance is the following: If a reliable estimate of the latter is available, then our theory 
essentially becomes parameter free. Furthermore, the closer the HF energy comes to the 
empirical binding energy of nuclear matter, the smaller the fluctuations are. 

II. BASIC FORMALISM 

While for the general presentation of the formalism we refer the reader to ref . [|I[] , here we 
confine ourselves to shortly revisit the basic equations we shall need in the following. These 
stem from the theory of nuclear reactions developed by one of us 0, where the projection 
operators P and Q are introduced and the Schroedinger equation is split into two coupled 
equations 



(E - H pp )(P^j) = H pq (Qi/j) 
(E-H QQ )(Q^)=H QP (P^), 

with obvious meaning of the symbols. From these the equation 
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for the system's wave function projected into the P-space follows. In (2.2) the standard 
definitions 
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have been adopted. 

The strategy of our approach is first to implement the above formalism by splitting the 
Hilbert space in a P sector embodying all the "gentle" physics and in a Q sector where the 
rapidly time dependent physics associated with the violent two-body short range collisions 
among nucleons is inserted. Clearly, an exact prescription for such a partition cannot be 
given: This, however, might be more an advantage than a flaw, since it adds to the flexibility 
of the method. 

Next, an energy averaging procedure is performed into the Q -space: Specifically, the 
fast energy-dependent operator 1/eg is averaged, leaving the weakly energy- dependent Wqq 
as it stands. 

Thus, the equation 
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is obtained. In (|2.4j ), Eq should be viewed as the MF energy and the angle brackets mean 
energy and phase averaging. Energy averaging is performed with the smoothing function 
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(2.5) 



e being a real, positive parameter, which moves the poles of the Q-space away from the 
P-space (the more so, the larger e is), thus quenching the influence of the former over the 
latter. 

With ( |2.5| ) one obtains, by applying the Cauchy theorem to a contour going along the 
real axis with a small semi-circle described positively around the E — e pole, 



j_\ j_ r dE x i 

v e Q / 2m J-oo E-(E -e)E- H QQ - W QQ {E) 
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E -e-H QQ - W QQ (E = E - e) " E Q - e - H QQ - W QQ (E = E ) ' 



(2.6) 



In performing the integration in the complex energy plane, account has been taken of the 
position of the poles associated with the operator 1/eg (they are located in the ImE < 
half-plane) ||. 

Furthermore, the approximation performed in ( [2.6|) exploits the mild energy dependence 
of Wqq, but of course it requires, to be valid, a not too large value of the parameter e. On 
the other hand, the latter cannot be too small in order to lessen the impact of the Q-space 
on the P one. Hence, e is constrained to vary in a rather restricted range. 

According to the above, we obtain for the P-space energy averaged Hamiltonian 7Y the 
expression 

ft = H PP + V pq Vq p ^ -, (2.7) 

the associated equation for the averaged wave function reading 

H(Pij) = E (Pij). (2.8) 



The energy En in fl2.8|) is commonly referred to as the MF ground state energy. In ( |2.7|) , in 
analogy with ref. M, the coupling operators 



V P Q — H P Q, 



En — e — E 
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K,i ~ ' ~ E -H QP (2.9b) 
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have been introduced. They should be viewed as representing what is left out of the inter- 
action among the nucleons after the average value of the latter has been removed. We note 
that V is zero when e vanishes. 

The above definitions are convenient because in terms of them the system fl2.ip can then 
be exactly recast as follows: 



As shown in [|TJ, is 



{E-H){Pi>) = V PQ {Qi>) (2.10a) 

(E - H QQ )(Q^) = V Q p(P^). (2.10b) 



E - E = (MVpQ v , V Q p\<h) = AE, (2.11) 

h> — lino 



with 

hQQ = H Q Q + V Q p f — _nj) VPQ 

= H QQ + W QQ . (2.12) 

In ( |2.12| ) the prime implies the dropping of the \<f>o) component of l/(E — H). 
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III. ENERGY FLUCTUATION 

The physical significance of ( |2.11| ) relates to the statistical fluctuations of the energy 
around the average value E (the MF energy). In fact, the average over the random phases 
of the wave-functions of the Q-space leads to 

(AE) = 0. (3.1) 

To implement the above constraint we split AE according to 

AE = J2 AR m, (3.2) 



m=l 



with 



and 



AE m = (^olVoiGiV^ • • • V m _i im G m V mim _i . . . V w \4> ), (3.3) 

(3.4) 

Gk = F—h v r v — (3 ' 5) 

& ~ n QQ ~ v k,k+l^k+l Vk+l,k 



V i:i = QiVQj. (3.6) 



Actually, in (|3.2|) the index m labels different sets of states of increasing complexity, which 
contribute independently to the error, the number of sets being r. Moreover, for each m 

(AE m ) = 0; (3.7) 

hence (|3.1|) follows. 

We shall confine ourselves to consider the first set (m = 1) only and will enforce ( |3.7| ) by 
defining 

AE! = (0ol VoidtWo) - [(0o|K)iG 1 V lo |0o)] A v , (3-8) 

whose square we shall evaluate in the following. 

Now, while (AE m ) vanishes, ((AE m ) 2 ) does not. We compute it, again assuming random 
matrix elements, as in |1[, but dropping a severe assumption there made. For this purpose 
let the ground state, either of nuclear matter or of a finite nucleus, be 0o- Assuming the 
residual interaction to be of two-body character, then the m — 1 set will be made up of 
2p-2h states ip2/3 and the propagator of the system in the m = 1 sector will obey the equation 

G^fof, = {E- h n )ikp = {E- e 2 ^2f3, (3.9) 

the index j3 embodying the quantum numbers identifying the 2p-2h states. Then 



(5E,) 2 = ((AE,) 2 ) - (AE 1 Y (3.10) 

_ /^ K^|V|0o)| 2 |(^ 7 |V|0o)| 2 \ /j, \(^\V\cf>o)\ 2 \ 2 
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+ E(^2/ 3 |^|0o)*^2 7 |^|0o))((^|^|0o)^2 7 |V r |0 O )*) 
/3 7 

+ fourth order terms — ( E I (V^l^l^o) | 2 ) 



where an average 2p-2h eigenvalue €2 has been introduced in the denominator, which has 
not been averaged owing to its smooth energy dependence. 

Next, neglecting in the above the fourth order and because of the randomness of the 
phases of the states, we finally arrive at the expression 

(^i) 2 = rF ^ 2 E(K^|V10o)| 2 ) 2 , (3.11) 

\& - e 2j g 



where the sum goes over the 2p-2h states contained in the energy range AE\ = e around e 2 
over which the average is carried out. 

IV. A SIMPLE PROJECTION OPERATOR 

In order to apply the formalism previously outlined a choice for the projection operator 
P should be made. We start with the simplest one, namely we take 

P=\Xhf){xbf\, (4- 1 ) 

Ixhf) being the HF ground state wave-function (the Fermi sphere in the nuclear matter 
case) . 

Then, by combining (|4.1|) with eq. (2.8), it is not difficult to derive the MF equation (see 



ref. [0) 

£0 = E HF + _ £ -, (4.2) 

£L/q — e — hi 

which relates the MF energy E , the HF energy E-^p and the true energy of the system E. 
In (|3) 
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(3 2 = £ |^ 2p - 2/i |V|xhf)| 2 . (4.3) 

2p-2h 

To proceed further, we elaborate the relation ( [3.11D in order to express the fluctuations of 
the ground state energy around E also in terms of /3 2 . For this purpose, we write 

v 2 ' /3 7 

l -— 2 A | J2 (<V^»I v|^,>*(v^|v|^,>) 1 E ((MvifaHMvifa)*) 



pi ) $i 



(E-e 2 

+^ I E {(^\V\<Po)(^\V\^r) \ E (^2/ 3 |V|0o)*^2 7 |^|0o)) 
1/37 J /?7 

= 2 .A I £ (<V^I V|^o>*<V^|V|^o>) 1 E ((^|^|0o)*(^2 7 |^|0o)). (4.4) 

^ ~ ^ I 07 J 07 

In the above, the operator „4 corresponds to the arithmetic mean of the quantity within the 
angle brackets over the set of the 2p-2h states in the energy interval e centered around e 2 - 
Moreover, because of the randomness of the wave functions, the sum over (3 ^ 7 vanishes. 
We thus obtain 

^|E((^|V|0o}*^2 7 |V|0o))| =^E(K^7Wo>| 2 ) 

= ^\ (4-5) 

M being the number of states over which the sum is performed. 

To obtain the latter we make use of Ericson's statistical theory of nuclear states ||, 
which yields for the density of the spin J nuclear levels set up with N particle (p)-hole (h) 
(iV=p+h) excitations the expression 

&& J) = i^7» expH2J + lM] ' (46) 

where 

is the so-called spin cut-off parameter. In (|4.7| ) a is a constant and T the system's moment 
of inertia. Note that the energy £ in ( |4.6|) is the positive excitation energy of the states of 
concern. 

The number of levels is then given by the expression 



A£l ^ (4) 



A/-=^ = ep^(£,0), (4.8) 



T>2 being the distance in energy between the spin-zero 2p-2h excitations, - the ones we 
confine ourselves, for the sake of simplicity, to consider. Hence, from ( |4.4|) , it follows 



((E - E f) 



P 4 
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and, by taking the square root, 



e 2 ) 2 AA 



(4.9) 



E — E 




(4.10) 



As in ref. [|IJ, eqs. Q4.2| ) and ( |4 . 1 0| ) together make up two non- linear systems, eq. ( |4 . 1 0| ) being 



taken with the plus and the minus sign, respectively, in the first and in the second system. 
The difference with the formalism of ref. [|1] lies in the factor 1/yM. 

In our approach the unknown quantities in the systems are either E and E Q , if a model 
for the matrix elements of the residual interaction is available, or Eq and (3 2 , if E is ex- 
perimentally measured and one wishes to extract the residual effective interaction from the 
data. In this paper, we shall follow the second option. For definiteness we call the solutions 
of the first system E l Q and E (or (3f) and, likewise, the solutions of the second system Eq 
and E u (or (3 2 ), the labels / and u referring to the lower and upper boundaries of the energy 
band embracing the fluctuations of the ground state energy of the system. 

Now, both systems lead to the quadratic equation 



El 



E (E HF + E + e) + E HF (E + e) - p 2 







(4.11) 



for the mean field energy. We take the solution with the minus sign in front of the square 
root, namely 



E = 1 {(£hf + E + e) - V /(£ HF -£-e) 2 + 4/3 2 } 



(4.12) 



since then, in the limit of vanishing residual interaction, which entails the vanishing of e (see 
eqs. Q2.9|)), we obtain 



En = E, 



(4.13) 



which is consistent with eq. ( 4.10 ) in the same limit. Then, when /3 2 — > the fluctuation of 
the ground state energy vanishes and the latter coincides with the mean field value. 

Concerning the residual effective interaction, one gets from the first system the expression 
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+ Ae(E l - E HF ) 



(4.14) 



and from the second system the expression 



0t = V-^/2 



-E u - ae 



E u {-1 + JN/2) - e(-l + JN/2a) + £ H f 



+ \ 



E u (-1 + t/jV/2) - e(-l + VJV/2a) + £ H f 



4e(£ u - E HF ) 



(4.15) 



In the above, in conformity to the notations of ref . |T) , we have set 

e 2 = ae 



(4.16) 



and the solutions with the minus sign in front of the square root have been discarded 
because they correspond to a residual effective interaction non-vanishing with e. Clearly, if 
an estimate of the average energy of the 2p-2h states can be derived, then a in (|4.16[ ) should 
not be viewed as a parameter, but merely as a convenient tool to express e 2 . 
Our solutions should fulfill two constraints: The first one is 



Ek 



E Q, 



(4.17) 



namely the two mean fields stemming from the two systems should be identical. The second 
relates to the spectroscopic factor S defined through the normalization of the ground state 
wave-function projected onto the P-space, namely 



{Pip\Pip) = S 2 



and found to be 
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(4.18) 



(4.19) 



in ref. QTJ . We require the values of S 2 , which approaches one (as it should) when (3 2 — > 0, 
to be "reasonable" when computed with our solutions for E and j3 2 on the upper and lower 
bounds of the energy band, respectively. Actually, as it is apparent from ( |4.19| ), S 2 depends 
also upon the energy derivative of the residual effective interaction. This reads 



M = 4^{ £ , (1 + V ^7i)_ e(1 + V ^ 
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[E^l + ^M/2) - e(l + yJN/2a) - Phf] 2 + 4e(E, - E w ) 
on the lower bound and 



(4.20) 



% = ^^U u (-1 + y/JT/2) - e(-l + y/JT/ia) + E UF 

+ J E u {-1 + yfif/2) - e(-l + 0\/y2a) + E RF 
+ (E U - at) 



+ Ae{E u - E nF ) 



(E u (-1 + JAf/2) - e(-l + jAT/2«) + E nF )UN/2 - 1) + 2e 



Af/2-1 

^[£„(-l + /A/"/2) - e(-l + ^jM/2a) + £ HF ] 2 + 4e(£ u - £ H f) J 

(4.21) 
on the upper one and reduces to the corresponding expression found in ref. [||] when J\f — > 2. 

V. NUCLEAR MATTER 

In this section we compute the nuclear matter mean field energy Eq together with its 
error. In addition we obtain, by solving the systems of eqs. (|4.2| ) and ( |4.10 ), the lower (f3f) 
and the upper ((3%) residual effective interaction. 

To carry out this task we exploit the empirical nuclear matter energy 

E NM = [-16 + 39.5(k F - 1.36) 2 ] MeV, (5.1) 

with kp in fm _1 . Notice that formula (|5.1|) embodies not only the experimental binding 
energy, but the saturation density and the compression modulus of nuclear matter as well. 
Furthermore, to start with we take E^ F as given in ref. |l|, where a simple, four parameter 
two-body force with short-range repulsion followed by intermediate-range attraction was 
found to yield a minimum for the HF energy of -7.3 MeV at kp = 1.78 fm _1 for some choice 
of the parameters. 

Concerning M we assume, for a nucleus of mass number A, the following values for the 
quantities entering into its definition: 

a = - MeV" 1 (5.2a) 

T = -m N I drr 2 p(r) = —m N p R 5 (5.2b) 

3 A , 

9 = ,—■ (5.2c) 

A 6p 

In the above, ( |5.2b|) corresponds to the moment of inertia of a constant mass distribution 
p(r) = po = 2kp/3n 2 of radius R = tqA 1 ^ 3 (r = 1.2 fm), m N being the nucleon's mass, and 
the Fermi energy e F is computed with the standard nuclear matter value of the Fermi wave 
number, i. e. k F = 1.36 fm -1 . 

Inserting the above in ( [4.6| ) we get, for a A = 208 nucleus taken as a good representative 
of nuclear matter, 

pfp-2h( E , J = Q) = 3.76 • 1(T 2 £ 9/4 exp(-0.0013/v / ^) MeV~\ (5.3) 
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£ = ae — -Enm being the excitation energy of the 2p-2h states. 

Next, by introducing the fluctuation energy W and setting accordingly 

W 

Ei = £ NM - — (5.4a) 



and 



W 

E u = £ NM + — ( 5 - 4b ) 



we obtain all the ingredients required to get, from (|4.14j ) and ( |4.15|) , (3f and /3 2 as a function 
of e and a for given values of W and £ . 



With these the equal mean field condition (|4.17|) can be cast into the form of the following 
equation 



E l - ^HF -Et- e)2 + 4(3? = E U - ^(E UF - E u - e) 2 + 4/3 2 , (5.5) 

which we have numerically solved for W — 1, 2 and 3 MeV and for £ — 14 MeV. The latter 
value entails e 2 = ae = -2 MeV and corresponds to £ = 2fnu, huj = Al/A 1/3 (= 7 MeV 
for A = 208) being the parameter of the harmonic oscillator well modelling the shell model 
binding potential. 

Our results are displayed in Fig. |l|a, where a continuum set of solutions is seen to exist 
for each value of W. Should the above estimate of £ be adequate, then the intercepts of the 
hyperbola a = — 2/e with the equal mean field curves of Fig. [i] would set the values (one 
for each choice of W) of the sole parameter entering in our theory, namely e, thus providing 
the consistent solutions of our approach. 

Of course, the estimate £ = 14 MeV is only orienting: Therefore, we display as well in 
Figs. [TJb and |I]c the curves corresponding to the choices £ = 18 (e 2 = 2) MeV and £ = 22 
(e 2 = 6) MeV. Note the dependence upon £ of the solutions of (|5.5|) . Indeed, the larger £ is, 
the smaller the energy fluctuation becomes, as it should be expected from the fluctuation 
equation ( |4.10| ). 

Finally, we collect our results, namely the values of e and a corresponding to the inter- 
cepts for W = 1, 2 and 3 MeV, together with the associated values of E l , Eq, S l , S u , f3f 
and Pl for £ = 14, 18 and 22 MeV in Tables |, and [In], respectively. 

These predictions of our theory should be examined with the guiding principle of having 
a parameter e neither too large, nor too small and a consistent spectroscopic factor S. 

From Tables [I], |TJ and [TO] we infer that the present theory comes indeed close to fulfilling 
the above requirements for a fluctuation energy, say, W ^ 1 -i- 2 MeV and an average 
excitation energy of the 2p-2h states £ < 20 MeV. 

Notably, this outcome is reached with a major reduction of the residual effective interac- 
tion with respect to the bare /3 2 , which, with the parameters characterizing the force taken 
as in |I[, reads (3 2 = 5.2 ■ 10 4 MeV 2 /nucleon. Thus the findings of ref. |I[ concerning both 
the size of the error of the nuclear matter ground state energy and the quenching of the bare 
interaction appear now established on a firmer ground. 

Concerning the spectroscopic factor, its value comes close to one, a result consistent 
with our effective interaction which is too weak to affect the occupation probability of the 
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FIG. 1. Solutions of eq. ( j5.5|) corresponding to a width of the energy fluctuations of 1, 2 and 3 
MeV (solid lines); the dashed lines correspond to ae = t2 = -Enm + £'■ (a) £ = 14 MeV; (b) £ = 18 
MeV; (c) £ = 22 MeV. 
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w 


e 


a 


rpl — rpu 

E = &o 


S 1 


S u 


A 2 


(MeV) 


(MeV) 




(MeV) 






(MeV 2 /nuc 


1 


2.56 


-0.78 


-16.11 


1.11 


0.91 


24.12 


2 


10.62 


-0.19 


-16.16 


1.06 


0.95 


109.35 


3 


23.42 


-0.085 


-16.25 


1.04 


0.96 


249.75 



Pi 

;V 2 /nuc 

35.25 

131.70 

283.56 



TABLE I. Solutions of eq. ( M ) and ae = -Enm + £ for £ = 14 MeV and for three values of 
the width W of the energy fluctuations; there are also reported the corresponding values of mean 
field, spectroscopic factors and effective interactions. The bare potential used in the calculations 
yields a HF field with a minimum of -7.3 MeV at kp = 1.78 fm _1 , a HF energy of -4.99 MeV at 
k F = 1.36 fm" 1 and /^ are = 5.2 • 10 4 MeV 2 /nucleon. 



w 

(MeV) 


€ 

(MeV) 


a 


rpl — rpu 

(MeV) 


S l 


S u 


A 2 

(MeV 2 /nucleon) 


2 

r J u 

(MeV 2 /nucleon) 


1 

2 
3 


8.15 
32.40 
71.49 


0.25 

0.062 

0.028 


-16.04 
-16.09 
-16.16 


1.04 
1.03 
1.03 


0.96 
0.97 
0.96 


85.11 
349.39 
783.36 


96.16 

371.58 
816.86 



TABLE II. As in Table |, but for £ = 18 MeV. 



w 

(MeV) 


e 

(MeV) 


El = E% 
(MeV) 


S l 


gu 


A 2 

(MeV 2 /nucleon) 


P 2 
(MeV 2 /nucleon) 


1 

2 
3 


19.78 

78.46 

173.85 


0.30 -16.02 
0.076 -16.05 
0.035 -16.12 


1.02 
1.03 
1.03 


0.98 
0.97 
0.96 


213.03 

858.10 
1919.00 


224.06 

880.24 

1952.38 






TABLE III. As in 


Table H 


[ but for £ 


= 22 MeV. 





W 


f 


a 


rpl — rpu 


S' 


S u 


(MeV) 


(MeV) 




(MeV) 






0.4 


3.99 


-0.50 


-16.01 


1.04 


0.95 


0.5 


6.24 


-0.32 


-16.01 


1.04 


0.95 


1 


24.80 


-0.081 


-16.03 


1.05 


0.94 



A 2 



ill 



(MeV 2 /nucleon) (MeV 2 /nucleon) 



14.19 
22.41 
91.10 



15.68 
24.28 
94.85 



TABLE IV. As in Table |, but with a bare potential yielding a HF field with a minimum of 
-13 MeV at k F = 1.5 fm" 1 , a HF energy of -12.28 MeV at k F = 1.36 fin -1 . 
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single-particle levels inside the Fermi sphere (in [[jj a sign error in d(3 2 jdE yielded quite small 
values for S). The discrepancy with current estimates, which yield S ~ 0.7 ||, is likely to 
stem from the restricted P-space we have adopted. Work is currently in progress to shed 
light on this point (by enlarging the P-space) |?J . 

Another important point, is whether our conclusions remain true for a different choice 
of Phf- Actually, the HF energy employed by us should rather be viewed as a tool needed 
in our framework than as a quantity precisely computed. Indeed, in order to calculate 
Phf? one should start from a realistic two-body force, deduce the associated G-matrix, thus 
accounting for the deterministic aspect of the short-range nucleon-nucleon correlations and 
then perform a HF (in fact, a Brueckner-Hartree-Fock) calculation. This we have not done. 
We have however varied the parameters of our simple two-body force to get Phf closer to 
the empirical energy, namely getting a HF minimum of -13 MeV at kp = 1.5 fm _1 with 
roughly the same empirical compressibility of nuclear matter, thus yielding Phf — —12.28 
MeV at k F = 1.36 fin -1 . 



The results thus obtained are quoted in Table IV. It is clearly apparent that the closer 



Phf is to the true value of the energy, the thinner the fluctuations are. 

VI. CONCLUSIONS 

In the present research we improve upon the statistical treatment of the mean field of 
the atomic nuclei and nuclear matter we have recently proposed |jj. Specifically, we deal 
more realistically with the equation expressing the fluctuation of the energy of the ground 
state by exploiting the Ericson's formula for the density of the 2p-2h states. The latter 
in turn allows us to introduce a self-consistency condition for the mean field energy, which 
correlates the value of the energy averaging parameter e with the average excitation energy 
S of the 2p-2h states. Should the latter be known, then no parameter would be left in our 
theory. 

Lacking such a knowledge, we have explored the dependence of our predictions for the 
ground state energy fluctuation W upon S. As expected, we have found that the larger S is, 
the more remote the Q-space becomes: hence, the smaller the W one gets. However, while 
£ on the one side cannot be smaller than 2fiu), on the other cannot be too large because 
its actual value relates to the average energy of those 2p-2h states connected to the ground 
state of H by the residual force V, which is found to be very weak. 

In principle also a large e would push away the Q-space, but then our theory, as previously 
shown, would not work any longer. In fact, as seen in the figures, the solutions associated 
with a large e also correspond to a large W, which does not make any sense. Rather, the 
trustable solution of our theory is the one corresponding to a W of one, at most two, MeV 
and to an energy averaging parameter e = 2 -S- 3 MeV, these outcomes being accompanied 
by a quenching of the bare interaction of more than an order of magnitude. 

The latter finding should not surprise: in fact, our residual interaction has been dove- 
tailed to vanish when e — > and e actually is about 2 -h 3 MeV. Since the physics associated 
with the Q-space extends up to several hundred MeV, the strong quenching of Vpq (Vqp) 
with respect to Hpq (Hqp) immediately follows. 

The above results, however, are obtained with a HF energy yielding a poor approximation 
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for the true energy. A E-^p closer to the latter entails, as we have seen, a smaller W (but 
not a smaller e). We argue accordingly that our theory provides an appreciation of the error 
of any theory attempting to calculate the ground state energy of nuclei and, of course, the 
nearer the theory comes to the true eigenvalue, the smaller the associated error is. 

Finally, the unrealistic, but consistent, value of our spectroscopic factor has already 
been commented upon. Here we wish to remark that the better estimate of S we expect by 
enlarging the P-space and by allowing configurations more complex than the 2p-2h ones to 
contribute, should go in parallel to a better appreciation of the error as well. 
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